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Abstract 



In a recent series of papers evidence has been presented for correlations between solar activity and nuclear decay rates. This 
includes an apparent correlation between Earth-Sun distance and data taken at Brookhaven National Laboratory (BNL), and at 
the Physikalisch-Technische Bundesanstalt (PTB). Although these correlations could arise from a direct interaction between the 
decaying nuclei and some particles or fields emanating from the Sun, they could also represent an "environmental" effect arising 
from a seasonal variation of the sensitivities of the BNL and PTB detectors due to changes in temperature, relative humidity, 
background radiation, etc. In this paper, we present a detailed analysis of the responses of the detectors actually used in the BNL 
0^ and PTB experiments, and show that sensitivities to seasonal variations in the respective detectors are likely too small to produce 
the observed fluctuations. 
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■1. Introduction 

' In a recent series of papers iH 0, S] evidence has been pre- 
sented suggesting a correlation between nuclear decay rates and 
solar activity. This evidence came from analyses of data from 
three independent sources. The first was data taken at Purdue 
University during a period which included the solar flare of 
2006 December 13, which exhibited a statistically significant 
'dip in the counting rate of ^"^Mn coincident in time with the 
solar flare fit]. The second was data from a measurement of 
the half-life (T1/2) of ^^Si at Brookhaven National Laboratory 
(BNL) and the third data set was an extended study of de- 
lector stability via measurements of the decay rates of '^^Eu, 
'^"^Eu, and ^^''Ra at the Physikalisch-Technische Bundesanstalt 
(PTB) in Germany The BNL and PTB data revealed pe- 

riodic variations in measured decay rates which approximately 
correlated with Earth-Sun distance. Taken together the data of 
Refs. d 0. S B B S] raise the possibility that nuclear decays 
are being influenced by the Sun through some as yet unknown 
mechanism, perhaps involving solar neutrinos. 

The ^"*Mn data presented in Ref. 01] were a series of succes- 
sive four-hour counts taken over the course of two months as 
part of a half-life measurement. The selection of a four-hour 
count time allowed for a time resolution capable of detecting 
the solar flare, which lasted approximately 43 minutes |7]. Due 
to the short counting period, the counts were presumably insen- 
sitive to typical seasonal environmental effects, such as temper- 
ature, pressure, background, etc. 
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However, the same cannot be said of the data of Refs. 
iSSSl since the annually varying decay rates reported in these 
references could in principle arise from a number of conven- 
tional systematic influences on the detectors employed in these 
experiments rather than from a modulation of the decay rates 
themselves. The object of the present paper is to analyze the 
BNL and PTB data in detail in an effort to disentangle envi- 
ronmental influences on the detector systems used in each ex- 
periment from a possible contribution originating from the Sun, 
which could affect the intrinsic decay rates. The detector sys- 
tem used in the ^"*Mn data set associated with the 2006 Decem- 
ber 13 solar flare will also be discussed. 

Although there are many mechanisms through which experi- 
mental detectors can be influenced by environmental factors, a 
useful starting point for the present analysis is the recent paper 
by Semkow, et al. j|8[] in which the authors attempt to explain 
both the BNL R "^Tand PTB ||1 0, El data in terms of sea- 
sonal temperature variations. In what follows we will analyze 
these experiments in detail and show explicitly why the mech- 
anisms proposed by Semkow et al. to explain the BNL/PTB 
data are not likely to be correct. This analysis can thus serve as 
a template for similar analyses of other past and future experi- 
ments on half-life determinations to evaluate possible seasonal 
environmental influences. 



2. The PTB Ionization Chamber Detector 

In the course of experiments at the PTB Q] which were in- 
tended to evaluate the stability of detector systems employed in 
the measurement of half-lives, e.g. '^^Eu (Ti/2=4936.6(20)d), 
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the authors reported an oscillation in the decay rate of the 
long-lived standard, ^^^Ra (Ti/2=1600(7)y) used in the mea- 
surements. These fluctuations were attributed by the authors 
to a "discharge eff'ect on the charge collecting capacitor, the 
cables and insulator to the ionization chamber electrode caused 
by background radioactivity such as radon and daughters which 
are known to show seasonal concentration changes L5J" . It was 
also proposed elsewhere [8] that the variations in the measured 
activity of the ^^''Ra source were caused by variations of the 
density of the argon working gas within the detector chamber 
itself, caused by seasonal changes in the ambient temperature. 
This and other potential systematic effects, such as changes in 
ambient pressure, humidity, background radiation or electronic 
fluctuations, are indeed possible explanations for the observed 
variations and should be examined. An alternative proposal to 
these explanations, as suggested in Ref. was that the os- 
cillations of the ^^^Ra measurements were not caused by local 
perturbations to the detector system, but were due to the influ- 
ence of fields or particles such as neutrinos emanating from the 
Sun. In what follows, we will examine possible environmen- 
tal influences on the detector system, and determine which of 
these, if any, could be responsible for the fluctuations in the 
^^^Ra measurements. 

To evaluate which environmental influences may have caused 
the fluctuations in the detector, it is helpful to understand 
how the detector utiUzed in the PTB experiments works. The 
measurements were taken with a high-pressure Any ionization 
chamber. Model IG12, manufactured by 20th Century (now 
Centronic, Ltd.), which was filled with argon to 20 atm (2 MPa) 
by the manufacturer, then welded shut prior to delivery to the 
user (Trevor Alloway, Centronic Ltd., private communication, 
2009). Therefore, the number of atoms of the working gas, ar- 
gon, within the detector is constant throughout the lifetime of 
the detector, irrespective of the external temperature, pressure 
or relative humidity. 

The ionization chamber for the PTB experiments was oper- 
ated in a "current" mode, as is typical with ionization chambers 
in this application. Here, the output current of the detector is 
generated by the collection of electron-ion pairs created by the 
ionization of the argon gas inside. These electron-ion pairs arise 
from photon interactions within the gas, and from the inter- 
action of photoelectrons and Compton electrons generated by 
photon interactions in the wall of the re-entrant tube, which then 
enter the working gas also causing ionization. Ionization cham- 
bers operate at sufficiently high bias voltages to prevent recom- 
bination of the electron-ion pairs generated by the interactions 
of photons, photoelectrons, and Compton electrons in the gas. 
In the case of the PTB measurements, that voltage was 500V. 
The electrodes (anode and cathode) of the detector are con- 
tained within the chamber of the detector itself, and these col- 
lect the electrons and ions created via the radiation interactions. 
The current output is then the result of the of electrons collected 
at the anode, which is ~ 1.44 X 10"^ Amp ■ Roentgen ' • hr ' 
(http://www.centronic.co.uk/ionisation.htm, 2009). One obvi- 
ous systematic variable that can be quickly examined is vari- 
ation in the bias voltage. However, ionization chambers are 
known for their insensitivity to applied bias voltage changes. 



since the ionization current is essentially independent of the ap- 
plied voltage in the ion chamber region 11 ill . This is due to 
the fact that, as pointed out above, the bias voltage applied in 
an ionization chamber is chosen to be sufficiently high to pre- 
vent appreciable recombination of the generated electron-ion 
pairs, while at the same time not providing enough energy to 
cause multiplication of electron-ion pairs within the detector It 
follows that the output current from the detector will be based 
solely on the number of electron-ion pairs generated by the in- 
cident radiation and collected at the electrodes within. 

Since it is the ionization of the argon gas within the detector 
that provides the measure of the energy deposited within by the 
incident radiation, and thus the activity of the sources being an- 
alyzed, it follows that changes in the properties of the gas within 
the detector, such as the atom density, could afffect the measure- 
ments. However, since the number of moles (n) of argon gas 
within the PTB detector is constant, the density n/V - p of the 
gas will also be constant, unless the volume (V) of the detector 
changes. Using the ideal gas law, PV - nRT, where R is the 
gas constant, T is the temperature, and P is the pressure, it fol- 
lows that if the volume is also constant (a non-constant volume 
will be addressed below) then the detector system is described 
by P\IT\ - P2IT2, where 1 and 2 denote the initial and final 
conditions, respectively. Hence, the only effect of a change in 
room temperature is to change the pressure inside the detector. 
And, since the ionization potential of argon is known to be a 
constant over all pressures from 1 to 200 bar |12], it follows 
that there should be no change in the response of the detector to 
changes in the room temperature if the volume of the detector 
is constant. 

We next consider the small effects that may arise from a vol- 
ume change due to expansion or contraction of the steel walls 
of the PTB detector as the room temperature increases or de- 
creases. The outer case and the bulkheads of the detector cham- 
ber are manufactured using Mild Steel (Trevor Alloway, Cen- 
tronic Ltd., private communication, 2009), and typical coeffi- 
cients of linear expansion, a, for steel range from -9x1 O^^K" ' 
to 17 X lO^^K"'. For our analysis here, we will use an average 
value of 13 X lO^^K"', and assume that the expansion of the 
metal is isotropic. Based on the stated dimensions of the de- 
tector, radius 9.25 cm and length 42.7 cm, and neglecting the 
small volume of the re-entrant tube (well), the initial volume, 
Vo, of the chamber at the initial temperature can be calculated 
to be ~11477.9 cm^. The change in volume, AV, arising from 
a change, AT", in temperature is 

VoH-AV = {I + aATf\{Li\ + aAT)] 

= y„(l+3aAr), (1) 

and hence the new volume of the cylinder will be 1 1478.3 cm^ 
for a 1°C increase in temperature in the room. The change in 
volume is thus ~4 parts in 100,000, and would then result in a 
fractional density change of -4 x 10"^, which is two orders of 
magnitude smaller than the density change proposed in Ref. Jsl] 
to explain the PTB data. Additionally, it should be noted that 
this is a decrease in density, not an increase as claimed in Ref. 
fs'l, since the same number of gas atoms now occupy a negligi- 
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bly larger volume. We can therefore conclude that temperature 
increases would lead to lower density of the working medium 
in the detector, and thus to lower counts in the summer, as ob- 
served in the PTB data. However, in order to achieve a 1 x 10"^ 
fractional change in density, which would cause a correspond- 
ing fractional change in the ionization and subsequent current, 
the room temperature would need to change by ~100°C, which 
is unrealistic. 

In contrast to the operation of the Centronic IG12 ioniza- 
tion chamber employed in Ref. ISJ and described above, whose 
density dependence is characterized p (T) w constant we next 
consider a hypothetical detector in which the pressure, P, and 
the volume, V, are held constant, which is the setup studied 
in Ref. lUt]. In this case, we find from the ideal gas law that 
pT » constant, so that Ap = -NT IT . The change in inten- 
sity A/ of y-rays passing through a thickness Ax of absorb- 
ing material with density p is given by A/// - jiphx where 
yu is the absorption coeflicient, and this leads to the famil- 
iar result [11] / = /oexp(-yupx). Without loss of generality, 
we can express the counting rate (C) as C = Co(l -I/Io), 
where Co is an appropriate normalization constant, and hence, 
AC/Ap = yuxexp (-yupx). To obtain the temperature dependence 
of C, we then replace |Ap| by |p| AT/T which gives 
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Although this reproduces Eq. 1 of Ref. fs*], we see from the 
preceding discussion that this result only applies to a detector 
for which pT - constant. This is not the case for the PTB 
detector, which is characterized by p{T) ^ constant as noted 
above. One can illustrate the practical difference between these 
two cases by using the values for yu, p, and x quoted in Ref. ^ : 
ji = 0.07703cmVg, p = 0.03567g/cm^ and x = 14.4cm. We 
find fipx = 0.0396, and hence, \AC/C\ = 0.98 |Ar/r| . By way 
of contrast, for the actual PTB detector, the only dependence 
of p(r) on T is through the small thermal expansion of the 
detector chamber, as discussed above. Using a coefficient of 
volume expansion y = 39 x 10 ''(°C) ' for steel, we find from 
the preceding analysis |Ap/p| = 0.98 \-AV/V\ = 4 x 10"^. For 
r=300K, we then find for the assumed temperature change in 
Ref. lH, AT = 0.91°C, \AC/C\ = 3.2 x 10"^ for the actual PTB 
detector. This is ~100 times smaller than the value |AC/C| = 
2.9 X 10"^ claimed in Ref. [8J from an inappropriate application 
of the ideal gas law to the actual PTB detector 

It follows from the preceding discussion that the expected ef- 
fects of temperature on the IG12 ionization chamber used in the 
PTB measurements are much smaller than the measured fluctu- 
ations in the ^^''Ra data. We then turn to other possible envi- 
ronmental effects, such as humidity, ambient air pressure, or 
ambient air density, which are less likely to affect the detector 
chamber itself since the detector is a closed system. Moreover, 
these effects are not likely to have an effect on the very small 
air gap between the source and wall of the detector well, since 
the Any ionization chamber is measuring the photons (gam- 
mas and bremsstrahlung) emitted by the sources, the fluxes of 
which would have little dependence on the marginal air density 
changes in the small air gap. 



We next turn to the possible contributions from background 
radiation of cosmic and terrestrial origins, such as the ^^^Rn 
variations considered by Siegert, et al. |5], and Schrader ifioll 
as the cause of the annual variation of the ^^^Ra measurements. 
Although ^^^Rn concentrations do exhibit seasonal as well as 
diurnal variations HI H [H [H [13, H [H, Q, the ac- 



tual activities are quite small, on the order of 15-90 Bq/m^, 
which are far below the activities of the sources being mea- 
sured. Specifically, an analysis performed by Abbady, et al. 
in Hannover, Germany [il8i1 (which is ~60 km from the PTB 
in Braunschweig), measured indoor concentrations of ^^^Rn 
inside the Centre for Radiation Protection and Radioecology 
(ZSR) at Hannover University. The data were collected for the 
period of one year, with measurements taken every two hours 
daily, from which monthly averages were calculated for each 
two hour interval. The data, shown in Table 1 of Ref. lll8ll . do 
exhibit a small annual variation in the concentration of ^^^Rn 
inside the ZSR offices (an institutional building which can be 
assumed to be somewhat similar to the institutional buildings 
at PTB), with the minimum and maximum concentrations oc- 
curring in February and May, respectively. Interestingly, the 
day/night variation shown in Ref. |18] appears to be larger 
than the annual variation, with the minimum and maximum 
occurring at 1600 and 0800, respectively. Plots of the time- 
averaged data and the monthly averaged data from Ref. |,18il 
are shown in Figs. [T]and|2] The diurnal eff'ect is clearly evi- 
dent, with the early morning maximum activity concentration 
possibly happening while the heating, ventilation and air con- 
ditioning (HVAC) systems are off or in a lower functional state, 
as is typical in large institutional buildings. It should also be 
noted that ^-^Rn concentrations are known to be dependent on 
air density, which is in turn dependent on temperature, relative 
humidity and pressure, as well as on precipitation. It follows 
that since the warmest part of the day is generally late after- 
noon, and the coolest is early morning, a diurnal effect is an 
understandable result. The annual variation is less obvious, but 
it appears that a maximum occurs in Summer/ Autumn. 

While the preceding information and analyses do not com- 
pletely rule out an effect caused by ^^^Rn, it appears that not 
only are the resulting activities quite small, but additionally the 
phase and period of the ^^^Rn concentrations do not match those 
of the ^^^Ra measurement variations reported in Ref. and 
later analyzed in Ref. [2J . 

Other contributions to background in the PTB measurements 
would be of terrestrial and cosmic origins, such as muons, or 
photons from '"'K decay. Efforts were made by Siegert, et al. 
to reduce the background by use of significant lead shielding 
around the detector (H , 21 1 . However, it should be noted that 
the ^^^Ra data provided by Schrader (H. Schrader, private com- 
munication, 2006) for the analysis presented in Ref. |Q] were 
already corrected for background, which obviates the need for 
analyses of these background effects. Nevertheless, a discus- 
sion of those possible backgrounds is presented below. An anal- 
ysis of environmental radiation at ground level was conducted 
by Wissmann at the PTB |19], utilizing data from the Cosmic 
Radiation Dosimetry Site, situated on a lake near the PTB labo- 
ratory. The results of the analysis, which corrected for fluctua- 
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Figure 1: 222Rn hourly averaged concentration obtained from Table 1 of Ab- 
bady et al. The data were taken at the Centre for Radiation Protection and 
Radioecology (ZSR) at Hannover University from July 2000 to June 2001. 

tions in weather and solar effects, exhibited an annual variation 
of +6.9nSv-h"\ which is equivalent to a dose variation equal to 
±1.92 X lO^^nSv ■ s ^ Given the initial amount of ^-''Ra in the 
PTB source, 300 /Jg, which would have had an initial activity 
of 11.1 MBq, the activity would be -88.8 MBq after 40 years 
due to the ingrowth of the radioactive daughters. Utilizing the 
more conservative activity value of 11 . 1 MBq, the gamma dose 
rate from the source can be estimated to be ~ 2. 1 x lO^^nSv • s"' , 
which is 7 orders of magnitude larger than the variation in the 
cosmic dose rate reported in Ref. [il9il . We can therefore con- 
clude that seasonal variations in cosmic radiation background 
would likely have little if any effect on the PTB detector sys- 
tem. 

Returning briefly to the background radiation from ^^^Rn, it 
should also be noted since ^^^Ra decays to ^^^Rn, the decay 
chain of ^^^Rn is identical to the decay chain of ^^''Ra. If the 
^^^Rn daughters are in equilibrium with the relatively constant 
^^^Rn concentration of x 55 Bq, the activity of ^^^Rn and its 
daughters at any time will be x 300Bq, which equates to a 
gamma dose rate of « 0.58nSv • s"' . This value is five orders of 
magnitude smaller than the dose rate of the ^^^Ra source itself. 
Hence, it can be safely assumed that even with a variation of 
the atmospheric ^^^Rn concentration of +100% (which is sig- 
nificantly larger than the actual variation noted in Ref. ifisll l 
the even larger activity of ^ 600Bq would yield a new gamma 
dose rate of ~ 1.16 nSv • s which is still four orders of mag- 
nitude smaller than the dose from the ^^^Ra source itself. This 
effectively rules out a discharge effect on cables, etc. arising 
from charged particles or photons produced by ^-^Rn daughters 
as the cause of the measured variations reported in Ref. |15|]. 

3. The PTB Solid State Detector 

The '^"^Eu and '^^Eu samples studied by Siegert, et al. fs*] 
were measured with three different semi-conductor detectors as 



Figure 2: Rn monthly averaged concentration obtained from Table 1 of Ab- 
bady et al. The data were taken from June 2000 through May 2001, and are 
arranged here January through December. See caption to Figure[T]for further 
details. 

part of the evaluation of detector system stability previously 
discussed. In the course of these measurements, a separate 
'^^Eu calibration source was used to measure the changes in 
efficiency of the Ge(Li) detector which was one of the three 
semi-conductor detectors examined. A plot of measurements 
of the efficiency of a Ge(L.i) detector versus time, which is pre- 
sented in Fig. 2 of Ref. ||5|], exhibits a similar period of 1 year, 
shifted by half a period (or about 6 months), with an ampli- 
tude of about 0.5%. The efficiency of a detector is determined 
by comparison of the measured counts taken with the detector 
to the expected output of a known source in a fixed geometry. 
Thus, for the Ge(Li) measurement of the '^^Eu point source, the 
minima of the efficiency curve represent a deficit in counts (with 
respect to the expected counts based on the known activity and 
half-life of the standard), while the maxima represent an excess 
of counts. Specifically, the minima of the efficiency plot fall 
near the beginning of the calendar year for the Ge(Li) detector, 
which is opposite to the signal seen in the ionization chamber 
measurements of ^^^Ra. This appears to call into question the 
conclusion of Jenkins, et al. that decay rates are simply de- 
pendent on the distance of the Earth from the Sun. However, 
the fact that the period of the '^^Eu data is the same, but shifted 
by half a period, still suggests a possible link to the Earth-Sun 
distance. 

This possibility is supported by an examination of the decay 
modes of '^^Eu which has two branches, K-capture (72.10%) 
and f3- decay (27.90%) \22L- The 1408 keV photon measured 
with the Ge(Li) detector Jsl] only arises from the K-capture 
branch, and is seen in 21.07% of the decays (there is no pho- 
ton of that energy in the decay branch) [22.1 . In contrast to 
the Ge(Li) detector, the ionization chamber is incapable of dis- 
criminating among the different photons measured, and hence 
could not track the 1408 keV photon separately, as it would be 
lost in the sum of all the photons emitted by the sample. The 
phase of the '^^Eu data could thus be understood if K-capture 
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modes responded to an external influence oppositely from fi^ 
decays. 

In Refs. ifllSl the conjecture was put forward that a field or 
particle emanating from the Sun, possibly neutrinos, might be 
enhancing or interfering with the decay of radioactive nuclides. 
Since the quantum mechanical details of the K-capture and (3^ 
decays are quite different, it is in fact possible that the response 
of the K-capture branch would be diff'erent from that of the yS~ 
decay branch. If these were in fact opposite effects, i.e. the de- 
cay rate for the decay branch would be increasing when the 
K-capture rate was decreasing, these competing eff'ects would 
tend to damp out any fluctuations. That such a possibility could 
actually occur is supported by data from ^^Mn (which also de- 
cays via K-capture) acquired during the solar flare of 2006 De- 
cember 13, flj- It was observed that the ^"^Mn count decreased 
from the expected rate during the three-day period encompass- 
ing the solar flare. Thus, if we assume that the solar neutrino 
flux increased during the solar flare, then the response of ^'^Mn 
and '^^Eu to increased solar flux would have the same phase 
in the sense that both decay rates decreased in response to an 
increase in solar neutrino flux. 

4. The BNL Proportional Counter 

We next turn to an examination of the environmental sensitiv- 
ity of the gas proportional counter and sample changing system 
used by Alburger, et al. f?] and Harbottle, et al. ^ to measure 
the half-life of ^^Si. The detector in question was cylindrical, 
~1.5 in. diameter, 2 in. depth, with a 1 in. end window, operat- 
ing on P-10 gas. The detector pressure was held constant by a 
device on the gas outlet vented to a barostated enclosure to limit 
fluctuations in the density of the P-10 gas. The detector itself 
was mounted on an automatic precision sample changer, which 
was entirely contained within the aforementioned barostated 
enclosure, as described in Refs. 141,19]. The window material 
was 0.006mm Kapton® (Dupont H-film), with gold vacuum- 
deposited on both sides (40-50 A on the outer surface, 120- 
150 A on the inner surface). A bias voltage of 2150 V was ap- 
plied to the detector, which was the voltage closest to the center 
of the beta counting plateau \A\. 

Gas proportional detectors utilize gas multiplication, where 
free electrons generated by the ionization of the gas by the in- 
coming particles are accelerated by the strong electric field (cre- 
ated by the bias voltage) to sufficient kinetic energy to cause 
further ionizations. Electrons freed in these ionizations can 
also be accelerated to cause additional ionizations, as long as 
their energy is greater than the ionization energy of the neu- 
tral gas molecules. This gas multiplication process, or cascade, 
is known as a Townsend avalanche, which terminates when all 



free electrons are collected at the anode 111 111 . With a well de- 
signed and maintained detector, the number of secondary ion- 
izations is proportional to the number of initial ionizations, but 
multiplied by a factor of "many thousands" ifTll . Almost all pri- 
mary ion pairs are formed outside the multiplying region, which 
is confined to a small volume around the anode, and the elec- 
trons from the primary ionizations drift to that region before 
multiplication takes place, "therefore, each electron undergoes 



the same multiplication process regardless of its position of for- 
mation, and the multiplication factor will be [the] same for all 
original ion pairs." 1 11 ] If the charge pulse of all of the electrons 
collected at the anode per event is larger than the discriminator 
level setting after amplification, a count is then recorded by the 
scaler 

The two sources used in the BNL experiment, ^^Si and "^^Cl, 
are described in detail in Alburger, et al.. Sec. 2 of Ref. 101 ■ The 
isotope ■'^Si undergoes y6~-decay (100%, Ep„,ax - 224.3 Ike V) 
to with a Ti/2=157(19) y p!^. The radioactive daugh- 
ter, ^-p (which is in secular equilibrium with ^^Si) undergoes 
y6--decay (100%, £/j™av^l710.5(21) keV) to (stable) with 
a ri/2=14.262(14) d lH. The isotope ^''Cl, on the other 
hand, decays via competing yS-decay branches (yS", 98.1(1)%, 
£/j,«fl.v=708.6(3)keV; /?+/e, 1.9(1)%) f2?]. Great care was taken 
by the BNL group in preparation of the ^'^Si source to ensure 
that no contaminants would be included that would confuse the 
data collection and analysis. Analyses of the spectra of each 
of the sources were conducted with a plastic scintillator to test 
for impurities or other problems, and the spectra were as ex- 
pected. The ^^Si betas were not noticeably evident in the spec- 
trum from the plastic scintillator, which was likely due to the 
effective thickness of the sample (17 mg/cm^) that served to de- 
grade and smear the low-energy betas, as noted by Alburger, 
et al. [4]. However, they did appear with the expected equal 
distributions when part of the same material was measured in a 
liquid scintillator The ^^Cl spectrum was as expected, 

as reported by Alburger, et al. [4] . 

The gas proportional detector described above was utilized 
in a differential counter system described in detail in Ref. 
which used a precision sample changing system that would 
alternate the sources during counting runs. The system was 
designed and built to allow maximum reproducibility, e.g. 
by utilizing precision micrometer heads which adjusted the 
source/detector distance to within +0. 001mm |9]. As stated 
before, this system was contained within a "pressure-regulated 
box in order to minimize the differential energy-loss effects that 
would occur with changes in barometric pressure" iJ]. Tests 
of the sensitivity of the counting system to changes in the box 
pressure, as well as to detector bias voltage, detector gas flow, 
and discriminator setting were carried out prior to the series of 
measurements, and the estimated effects on the ^^Si counts and 
the Si/Cl ratio are detailed in Table 1 of Ref. [4]. 

Notwithstanding the tests described above, Alburger, et al. 
identified an unexpected periodic fluctuation in the Si/Cl ratio 
with a period of approximately 1 year. Having virtually elim- 
inated other systematic variables, or at least having controlled 
them to the extent possible, Alburger's BNL group was left with 
temperature and humidity as possible explanations, which are 
known to change the density of air (even if the pressure of the 
box containing the detector system is held constant). Increases 
in temperature will decrease the density of air, as will increases 
in relative humidity, although humidity has a much smaller ef- 
fect than temperature. To summarize, warmer more humid air, 
i.e. summer air, is less dense than cooler, drier winter air. This 
is depicted in Fig. |3] which plots the temperature and relative 
humidity dependence of the density of air, based on equations 
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developed and recommended by the Comite International des 
Poids etMesures (CIPM) ||27l |28l|29ll . 




Temp (°C) 



Relative Humidity (%) 



Figure 3: Plot of air density as a function of temperature and relative humidity. 
The data were calculated based on the formulae presented in Refs. | 27, 28, 231 



suggested by Alburger, et al. for temperature (T) (73+5°F) and 
relative humidity (RH) (35-76%), and assume that the maxi- 
mum density would be associated with the minimum T and RH 
(68°F, RH=35%), and the minimum air density with the max- 
imum T and RH (78°F, RH=76%). The air densities are then 
found to be l.lSSmg/cm-' and 1.155 mg/crn^, respectively. The 
energy-dependent linear ranges of energetic electrons for each 
condition set are shown in Fig. |4] which gives the linear range 
of betas of various low energies in air We have neglected the 
effects of the energy losses due to transport through the detec- 
tor window or sources, which are assumed to be independent of 
temperature. It should be noted again, that the ^^Cl source was 
4mm from the detector window, and the ■'^Si source 1mm away. 
Fig. |4]confirms the conclusions of the BNL group that changes 
in T and/or RH would have been too small to explain the fluc- 
tuations in their data. These conclusions are further supported 
by an MCNP analysis which we have carried out, the results of 
which will be described below. Before turning to this analysis, 
we first consider the remaining question of backgrounds which 
could impact the BNL detector system. 



The Alburger group at BNL did not begin collecting data on 
temperature and humidity until the last five months of the four 
year experiment, and during that time saw a range of tempera- 
tures of 72.4-74.7°F, and a relative humidity range of 35-76% 
10] • However, further investigations into data collected in other 
laboratories within the same building suggested that the aver- 
age temperatures remained well within the range of 70-76°F 
year round. The BNL group noted that the fluctuations appear 
to follow the same annual cycle as outdoor temperature. They 
also noted that if the indoor temperature and relative humidity 
were to track the outdoor conditions, then the data appeared to 
indicate that the fluctuations were environmentally based. The 
cooler, drier, and denser winter air could attenuate the lower 
energy ^*C1 betas, and would thus affect the ratio Si/Cl because 
the lower energy ^^Cl counts would be more weather depen- 
dent. However, Alburger, el al. concluded that "in order to 
produce the variations of ±3 standard deviations, the large hu- 
midity changes would have to be combined with temperature 
variations over a range of at least +5°F, which is larger than 
the probable actual range. We therefore conclude that system- 
atic periodic variations are present but that they cannot be fully 
accounted for by our tests or estimates" [4]. 

The conclusions of the BNL group can be supported by es- 
timating the effect of temperature and relative humidity on the 
range of betas over the spectrum of energies that would be seen 
from both sources. The range of an electron can be estimated by 
utilizing an analysis by Katz and Penfold 1^30^, which provided 
a range-energy (R - £„) relationship given by: 

R (mg/cm^) = 412£l-256-0.09541n£„ 

The linear range can be derived by dividing R by the density 
of the propagating medium, air, and expressing the electron en- 
ergy, £0 in MeV (it should be noted that these linear ranges are 
not necessarily straight line paths). We utilize the variations 




f," Energy (MeV) 

Figure 4: Linear range of betas in air as a function of beta energy. The graphs 
are calculated using the results of Katz and Penfold in Ref. 1 30]. The air densi- 
ties shown are based on the maximum and minimum values of temperature and 
humidity quoted by Alburger, el al. in Ref. Ql 

The BNL group measured an initial background counting 
rate of 6.6(3)/min. This background measurement was repeated 
after the conclusion of all runs by counting with the ^^Si source 
removed from its holder, and the '"'CI still in place (to ensure 
that there was no cross-talk between the sources). The rate was 
found to be 6.5(3)/min. Since the average count rates for the 
measurements were initially substantially larger, 21,500/min 
and 14,800/min for the ^^''Cl source and ■'^Si source, respec- 
tively, the background count rates were considered insignifi- 
cant, and hence were ignored by the BNL group. Additionally, 
since the sources were alternated during the counting procedure 
for each run, 30 minutes each for 20 cycles, the background 
would likely be the same for both of the isotopes. 
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5. Monte Carlo Modeling of the BNL Detector 

The Monte Carlo radiation transport software package, MC- 
NPX 113 in . was used to determine the effects of changes in 
ambient air temperature on the proportional detector system 
used by Alburger et al. f?]. The air density between the 
source and detector window was varied over the temperature 
range 40°F-90°F by use of the ideal gas law, corresponding 
to an approximate density range between 1.1350xlO"-'g/cm-' 
and 1.2486xl0"^g/cm^ for dry air. While humidity also af- 
fects total air density, this effect is presumed to be small (see 
Fig. [3]) relative to the effects of temperature change, and hence 
was ignored in this study. Additionally, the 40-90°F range that 
was analyzed encompasses both the minimum and maximum 
air densities expected for a wide range of both temperature and 
humidity in the laboratory. The physical parameters of the de- 
tector system (including the proportional detector, ■'^Si source, 
and ^^Cl source) were obtained from the literature f?, the 
BNL group's experimental notes (private communication, D. 
Alburger and G. Harbottle, 2005), and measurements of both 
the actual sources used in the experiment and a proportional 
detector of similar design. 

The cylindrical end-window proportional detector was mod- 
eled as a stainless steel (type 304) cylindrical shell with inner 
diameter 3.96875 cm, length 5.08 cm, and 3 mm wall thickness. 
The detector window was modeled following Harbottle et al. 
1^ as 0.006 mm Kapton® plated with vacuum-deposited gold 
on both the inner (125 A) and outer (45 A) surfaces. The detec- 
tor volume was filled with P-10 gas (90% argon, 10% methane 
by volume) with a calculated density of 1.6773xl0"^g/cm-'. 
This gas density within the detector was held constant. 

Due to the possibility that the lower energy betas from ^''Cl 
(708.6keV) were more easily affected by changes in air density 
than the higher energy ^^Si-^^P betas (1710. 5keV), each source 
was modeled separately, which also allowed us to account for 
the difference in proximity of each source to the detector. The 
■'^Si source, which was modeled in detail as described by Al- 
burger et al. 101, consisted of a brass base 3.18 cm in diameter 
and approximately 0.3175 cm (1/8 in.) thick with a 1.91 -cm- 
diameter, 0.8-mm-deep recess. The recess contained 47.7 mg 
of SiOa covered by 9.263x10"'' cm of aluminum foil. The MC- 
NPX electron source volume was defined uniformly over the 
volume occupied by the SiOa. According to Alburger et al., 
the ^^Si02 source was created approximately 13 years prior to 
the start of the half-life experiment [4], and hence ^^Si was cer- 
tainly in secular equilibrium with its radioactive daughter, ^^P, 
during the experimental time frame. The source energy spec- 
trum used in these MCNPX models consisted of equal contri- 
butions from the ^"Si (E,„„,. = 225 keV) and ^^^p (£,„„., = 1710.5 
keV) beta energy spectra calculated according to Fermi theory, 
which allowed us to incorporate the loss of the much lower en- 
ergy betas from ^'^Si. 

As was the case for the -^^Si source, the ^^Cl source was 
placed on a brass base 3.18 cm in diameter and 0. 3 1 75 cm thick, 
but without a recess machined into the top. The ^^Cl was spread 
over a diameter of 1 .9 1 cm in the center of the top surface of the 
base. Thus, the MCNPX electron source was defined as a flat 



disk 1.91 -cm in diameter on the surface of the brass base with 
an energy spectrum equal to that of ^^Cl {E,„ax - 709.2 keV) 
calculated according to Fermi theory. 

According to the BNL experimental notes, the distance of 
each source from the detector window was adjusted in order 
to obtain similar count rates from each, and this was achieved 
by having the ^^Si source 1.000 mm from the detector window, 
while the ^^Cl was placed at a distance of 4.000 mm. The source 
position in the MCNPX models for each source reflected these 
distances. 

As stated before, Alburger, et al. reported that the tem- 
perature variation was well within the range 70-76°F. We found 
that tally values varied so little over this range that they were 
statistically constant to within the minimum achievable relative 
error in MCNPX. Thus, a large temperature range of 40°F to 
90°F was chosen to improve the statistical determination of the 
function relating each tally to air temperature (and associated 
density). The results presented here were obtained using 1x10^ 
(for ^^Cl) or 2 X 10^ (for ^^Si-^^P) source electrons for each den- 
sity variation, and this provided relative errors of ~ 10 '* for the 
particle/energy current tallies, and ~ 2 x 10"^ for the energy 
deposition tallies. 

A number of MCNPX tallies were used to test the assertion 
that variations in air density with temperature would result in an 
associated increase or decrease in the probability that an elec- 
tron with a given energy will traverse the air gap between the 
source and detector window. A linear least-squares fit was ap- 
plied to the results of each tally as a function of air temperature, 
and these results are summarized in Table 1 . 

Of particular interest are the temperature dependent varia- 
tions in the energy deposition within the detector (Fig. |5]l , the 
electron current at the detector window (Fig. |6]l , and the en- 
ergy current at the detector window (Fig. |7]i , all of which can 
be used as surrogates for detector count rate. For each of these 
tallies, ratios of the value per source electron for ^^Si to that of 
^^Cl were calculated. Energy deposition, electron current, and 
energy current ratios were found to vary by 1 .77(40) x 10"^/°F, 
1.46(21) X 10-V°F, and 3.52(48) x 10"^°?, respectively. These 
results are echoed by the tally variations as a function of tem- 
perature for ^^Si-^^P and ^^Cl individually in Table 1. As ex- 
pected, all individual taUy values increase with increasing tem- 
perature due to the associated decrease in air density between 
the source and detector In all cases, the rate of change is greater 
for ^^Cl due to its lower energy spectrum (relative to ^^P) and 
greater source-detector distance. In general, detector events 
arising from the ^^Si source are due almost entirely to ^-P be- 
tas because the low-energy ^^Si betas are largely attenuated by 
the source itself I'T]. This explains the negative slope associated 
with the Si/Cl ratios as temperature increases. 

The ratios of electron current and energy current across a 
plane just above the source, and above the detector window 
were also examined (see Figures [8] and |9] ). These ratios di- 
rectly address the assertion that changes in air density result in 
count rate variations due to increased or decreased absorption 
of electrons in the air gap between source and detector The ^^Si 
electron and energy current ratio dependence was found to be 
6.8(32) X 10-^/°F and 5.1(33) x IQ-^TF, respectively, and for 
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Table 1 : Table of MCNPX results showing the sensitivities of the system per degree Fahrenheit. 





^^^Si- 
Per Source 


32p 

Norm.(70°F) 


H6 

Per Source g- 


CI 

Norm.(70°F) 


Ratio 


Det. E Deposition 
(MeV/ptcl) (A/°F) 


0.051(24) X 10-'' 


10.1(67) X 10-" 


0.34(40) X 10-" 


39.8(67) X 10-" 


-17.7(40) X 10-" 


Det. Window Current 

(A/°F) 


4.39(86) X 10-" 


12.1(33) X 10"" 


20.3(13) X 10-" 


35.3(33) X 10-" 


-14.6(21) X 10-" 


Det. Window E Cur- 
rent(MeV/ptcl) (A/°F) 


1.48(45) X lO-*" 


7.8(33) X 10-" 


4.32(32) X 10-" 


32.4(34) X 10-" 


-35.2(48) X 10-" 


Det./Source Current 
Ratio (A/°F) 


6.8(32) X 10-" 


7.1(48) X 10-" 


29.5(24) X 10-" 


41.5(48) X 10-" 


-46.5(64) X 10-" 


Det./Source E Current 
Ratio (A/°F) 


5.1(33) X 10-" 


5.2(47) X 10-" 


24.9(25) X 10-" 


33.7(48) X 10-" 


-37.9(63) X 10-" 
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Figure 5: The ratio of ^^Sito ^'^Cl energy deposition in the proportional detec- 
tor. Dashed lines indicate the 95% confidence interval for the calculated linear 
least-squares fit of the data. 



Figure 6: The ratio of -^^Si to ''Cl electron current (number of electrons) across 
the proportional detector window (into the detector). Dashed lines indicate the 
95% confidence interval for the calculated linear least-squares fit of the data. 



^"Cl, 2.95(24) X 10-V°F and 2.49(25) x 10-V°F, respectively. 
These temperature dependencies are all approximately 2 orders 
of magnitude smaller than the periodic variations observed by 
Alburger et al. 

As we noted previously, following a number of experimental 
tests, Alburger et al. concluded that it was possible that part 
of the 0.3% annual variation in observed count rates could be 
attributed to changes in ambient temperature (between 70 and 
76 °F), and relative humidity (35 to 76%) [4]. However, the re- 
sults obtained in this Monte Carlo analysis indicate that a tem- 
perature variation much larger than 6°F would be required to 
generate the observed count rate oscillations. Given that the ra- 
tio of energy deposition in the detector varies by approximately 
1.2 X 10-^/°F, the temperature of the air separating the source 
from the detector would have had to vary by well over 100 °F 
to account for the observed data. We therefore conclude that 
changes in temperature and relative humidity alone could not 
have caused the observed periodic count rate variations. 



6. Further Analysis 

In Fig. [To] , we exhibit the count rate ratio Si/Cl = 
N (Si) /N (CI) as a function of time, along with the outdoor tem- 
perature variation for the area surrounding the BNL |32] and 
1 /R^. We see that Si/Cl approximately correlates with 1 /R^ and 
anti-correlates with temperature, in both cases with observable 
phase shifts. Since at this stage we have no theoretical reason to 
suppose that Si/Cl should correlate (rather than anti-correlate) 
with 1 /R^, both of the correlations in Fig. [10] are significant. 
However, we have already ruled out temperature as the primary 
explanation of the Si and CI individual fluctuations, and hence 
we wish to understand whether influences emanating from the 
Sun, such as neutrinos, could be compatible with these data. 

In such a picture, the effect of analyzing the ratio Si/CI is to 
largely cancel out small fluctuations in voltage, currents, tem- 
perature, background, etc., which are common to both A^(Si) 
and (CI), while leaving intact the presumably larger effects 
due to the variation of 1/R^. From this perspective, the fact that 
the Si/Cl ratio does not exactly track 1/R^ is understandable in 
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Figure 7: The ratio of ^^Si to '*C1 energy current (MeV) across tlie proportional 
detector window (into the detector). Dashed lines indicate the 95% confidence 
interval for the calculated linear least-squares fit of the data. 



terms of the differing sensitivities of Si and CI to external per- 
turbations: the same nuclear structure mechanisms which en- 
dow Si and CI with very different lifetimes (172y and 301,000y, 
respectively) could presumably result in different responses of 
these nuclei to small perturbations. Evidence of this effect was 
first pointed out by Ellis |33], who saw an annual periodic sig- 
nal in ^^Mn, but not in '^^Cs when both nuclides were mea- 
sured for several years on the same detector system. Therefore, 
a long-term periodic signal due to the variation of l/R^ would 
survive in the ratio of Si/Cl. This can be understood quantita- 
tively by writing 



A^(Si) 

A^(Cl) 



A(Si)[l +e(Si)coswf] 
A(C1)[1 -He(Cl)coswf] 

^ ^^^^ ll + [e (Si) - e (CI)] cos cot] . 



A (CI) 

Here A, B, and e are constants, with e(Si) e(Cl), and e <«c 1 
has been assumed. In summary. Fig. [TO]is compatible with a 
picture in which any differences between Si/Cl and l/R^ are due 
to the varying sensitivities of Si and CI to an external perturba- 
tion coming from the Sun, as seen in Ellis' data. Although a dif- 
ferential temperature effect could in principle also explain Fig. 
[TOl this would require temperature-dependent effects which are 
much larger than allowed by the preceding analysis. 

7. Gravitational and Related Influences 

The preceding discussion leads naturally to the question of 
whether the observed effects in the BNL and PTE data could 
arise from a change in the gravitational potential of the Sun at 
the respective detectors as the Earth passes from perihelion to 
aphelion and back. In the framework of conventional General 
Relativity (GR) all clocks at a given space-time point run at the 
same rate, irrespective of their internal mechanism. In terms of 
Eq. |4]this would correspond to setting e(Si) = e(Cl), in which 



Figure 8: The ratios of electron current exiting the source toward the detector 
and electron current across the proportional detector window (into the detector) 
for -^^Si and ^'Cl. ^'Cl data have been shifted up by 0.2465 to facilitate visu- 
alization of both curves. Dashed and dotted lines indicate the 95% confidence 
interval for the calculated linear least-squares fit of the data for ""^Si and ^^Cl, 
respectively. 



case there would be no fluctuation in the ratio Si/Cl. Since this 
is in obvious disagreement with the BNL data, conventional GR 
can be ruled out as an explanation of the observed effects. 

However, this conclusion does not necessarily hold for vari- 
ous alternatives to GR. Will jjill has shown that in non-metric 
theories, clocks of different construction can in fact behave dif- 
ferently in a gravitational field. Although such a mechanism is 
possible in principle, it fails in practice on quantitative grounds. 
The Sun's gravitational potential O = GMq/Rc^ at the Earth 
is <£ = 9 X 10"^, and the fractional change AO between peri- 
helion and aphelion is AO = 3 x 10"^°. An effect this small 
would be undetectable given the BNL and PTB statistics, even 
if the coefficients multiplying AO in a non-metric theory were 
relatively large. Additionally, constraints on parameters which 
measure deviations from GR are sufficiently restrictive [i35ll to 
preclude the possibility of explaining the BNL and PTB data in 
terms of any known alternative to GR. 

A class of gravity-related theories which cannot be excluded 
at present are those in which scalar fields are introduced to in- 
duce a time variation in fundamental constants such as the fine 
structure constant a - e^/fic. We have discussed elsewhere 
llsll the possibility that although the simplest such theories can- 
not account for the BNL or PTB data, theories in which two or 
more scalar fields are introduced which influence both a and the 
electron/proton mass ratio (and possibly other quantities) might 
work. 



8. Environmental Influences During the Solar Flare of 2006 
December 13 

We turn in this section to a discussion of possible environ- 
mental influences on the Nal detector system used in y| during 
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Figure 9: The ratios of energy current exiting tlie source toward tlie detector 
and energy current across the proportional detector window (into the detector) 
for ^-Si and ^'Cl. The ^^Cl data have been shifted up by 0.2390 to facilitate 
visualization of both curves. Dashed and dotted lines indicate the 95% confi- 
dence interval for the calculated linear least-squares fit of the data for '"Si and 
^*C1, respectively. 



the 2006 December 13 solar flare. This is motivated by the po- 
tential connections between the ^'*Mn decay data collected in 
the time period surrounding the flare, and the BNL and PTB 
decay data which are the primary focus of the present paper. As 
we have already noted in the Introduction, the short duration of 
the solar flare allows us to rule out seasonal temperature, hu- 
midity and pressure influences on the detector system as poten- 
tial explanations of the observed decrease in the ^^Mn counting 
rate. Hence the focus of this section will be on other possible 
"environmental" influences on the detector system during the 
solar flare, such as neutrons, muons, etc. 

Although cosmic ray particles are known to be capable of 
influencing some detectors, there is a generic argument against 
any mechanism which attributes the drop in the ^'^Mn count rate 
to cosmic rays. This is based on the observation that the signal 
for a ^'*Mn decay in our system is the detection of the charac- 
teristic 834.8 keV photon resulting from the electron-capture 
process. 



e- + gMn(7'' = 3+) ^ v, + ^^Cr(7^ = 2^) 

^ llCr(g.s.)+y(834.8keY) (5) 



Hence the 834.8 keV photon which signals the electron-capture 
process is uniquely associated with the specific transition of the 
excited from the excited j'' = 2^ state to the 7^ - 0^ 
ground state (g.s.). The data points utilized in Ref. [1 ] represent 
the integral count within the Region of Interest (ROl) set on the 
834.8 keV photon for a four hour live-time. Even if there had 
been a large change in the cosmic ray flux of neutrons or muons 
during the flare, the chance that the nuclear or electromagnetic 
interactions induced by these particles would have accidentally 




1/1/82 7/2/82 1/1/S3 7/2/83 1/1/84 7/2/84 12/31/84 7/2/85 1/1/86 7/2/86 
Date 

— Si/C!3ptavg --USNOl/R'^2 — Outside T 

Figure 1 0: The counting rate ratio N (Si) /N(Cl) for the BNL data as a function 
of time, along with 1 /R^ and the outdoor temperature measured by NOAA near 
BNL. The correlation with l/R- and the anti-correlation with temperature are 
both evident. See text for further discussion. 



produced a photon at 834.8 keV, or deposited the energy equiv- 
alent of 834.8 keV within the Nal crystal, is extremely small. 
Moreover, the fact that a deficit of counts appeared in the ROI, 
and not an increase, indicates that this already very remote pos- 
sibility is highly unlikely. It is also unlikely that counts were 
selectively removed from the ROI as a result of changes to the 
cosmic ray flux. An alternative explanation might be that an 
increase in the cosmic ray flux would generate a large distribu- 
tion of energy depositions across the energy spectrum, and thus 
change the system dead time, leading to a loss of counts within 
the ROI. Dead time was tracked as a variable, however, and did 
not show changes over the course of the flare when compared to 
the counts preceding or following the flare period. Additionally, 
inspection of the spectra files also show no significant changes 
in the counts of channels above or below the 834.8 keV ^'*Mn 
peak. 

The above argument is bolstered by observations made by 
other groups during the 2006 December 13 flare. Although 
an analysis of data from world-wide neutron monitors by 
Butikofer, et al. and Plainaki, et al. SIstIi exhibited a -90% 
increase in the 1 -minute data at Oulu and Apatity, the overall 
neutron flux was still so smaU as to have a negligible effect on 
the detector system. This is particularly true in light of the fact 
that the Nal detector is insensitive to neutrons, and the flux of 
neutrons was too low to affect the ^'^Mn sample directly. 

Similar arguments allow us to exclude a change in the flux of 
cosmic ray muons resulting from a Forbush decrease |38] as a 
possible explanation of the flare data. A Forbush decrease is a 
rapid change in the flux of cosmic rays resulting from plasma 
clouds emitted by the Sun during a solar storm [381 . Evidence 
for a Forbush decrease during the 2006 December 13 flare ob- 
tained from the muon telescopes at the Basic Environmental 
Observatory (BEO) in Bulgaria and from Nagoya, Japan are 
shown in Figs. 5 and 6 of Ref. [38], respectively. The BEO data 
exhibited a sharp Forbush decrease beginning at ~15:30-16:30 



10 



UT on 2006 December 14, and reached a maximum at 05:00 
UT on December 15. The beginning of the Forbush decrease 
thus occurred 37 hours after the dip in the ^"^Mn counting rate 
reported in Ref. [JJ, which coincided in time with the solar flare 
at ~02:40 UT on December 13. Additionally, the minimum in 
the muon count rate associated with the Forbush decrease was 
observed more than 2 days later The data from the BEO are 
supported by data obtained at Nagoya, which exhibited a simi- 
lar time-dependence |38]. Given the timing of the Forbush de- 
crease relative to the dip in the ^"^Mn counting rate we can thus 
reasonably conclude that the dip was not likely the result of the 
response of the detector system to the Forbush decrease. 

The preceding arguments can be further strengthened by not- 
ing from Ref. jlj] that the ^"^Mn count rate started to decrease 
~1.7 days before the flare. It follows that the beginning of the 
decrease in the ^'^Mn count rate preceded the beginning of the 
Forbush decrease by more than 70 hours, and hence was not 
caused by it. 

In summary, we can rule out environmental explanations of 
the correlation between the solar flare of 2006 December 13 
and the observed dip in our ^'^Mn data |1] on several grounds. 
These include, (a) the observation that the solar flare lasted ap- 
proximately 43 minutes, which is too short a time for seasonal 
environmental influence to have affected the detector system, 
(b) the failure of these mechanisms to produce the character- 
istic 834.8 keV photon arising from ^"^Mn electron capture, or 
to account for the ~1.7 day precursor signal observed in the 
^"^Mn decay data; (c) the timing of the Forbush decrease which 
followed, rather than preceded, the dip in the ^'^Mn data. Addi- 
tionally, the observed changes in cosmic ray fluxes would have 
been too small quantitatively to account for the ^"^Mn data. 

9. Discussion and Conclusions 

As we noted in the Introduction, there are at present two 
competing general explanations for the apparent fluctuations 
observed in the BNL and PTB data: (a) they arise from the 
responses of the respective detector systems to seasonal varia- 
tions in temperature, pressure and humidity, and possibly other 
factors such as radon buildup, (b) the fluctuations arise from 
the decay process itself, due to some as yet unknown influence 
possibly originating from the Sun. By modeling the respective 
detector systems in detail, we have shown here that alternative 
(a) is an unlikely explanation for the observed BNL and PTB 
data, and hence by implication alternative (b) is more likely to 
be correct. The resulting inference, that nuclear decay rates are 
being directly influenced by solar activity is further supported 
by our analysis of ^"^Mn decay data acquired during the solar 
flare of 2006 December 13, which we have also analyzed. Our 
conclusion that the dip in the ^"^Mn data, which was coincident 
in time with the flare, was not likely attributable to changes in 
the cosmic ray flux during the flare, further strengthens the case 
that nuclear decays are being directly affected by solar activity. 
This inference can evidently be checked in a variety of experi- 
ments such as those described on Ref. Jstl, some of which are 
already in progress. 
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